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ABSTRACT 

Dynamic velocity dispersion and mass estimates are given for a sample of five X-ray luminous rich clusters of 
galaxies at intermediate redshifts (z ~ 0.3) drawn from a sample of 39 clusters for which we have obtained gravi- 
tational lens mass estimates. The velocity dispersions are determined from between 9 and 20 redshifts measured 
with the LDSS spectrograph of the William Herschel Telescope, and virial radii are determined from imaging 
using the UH8K mosaic CCD camera on the University of Hawaii 2.24m telescope. 

Including clusters with velocity dispersions taken from the literature, we have velocity dispersion estimates for 
12 clusters in our gravitational lensing sample. For this sample we compare the dynamical velocity dispersion 
estimates with our estimates of the velocity dispersions made from gravitational lensing by fitting a singular 
isothermal sphere profile to the observed tangential weak lensing distortion as a function of radius. In all but two 
clusters, we find a good agreement between the velocity dispersion estimates based on spectroscopy and on weak 
lensing. 

Subject headings: cosmology: observations — dark matter — galaxies: clusters: individual (Abell 914, 

Abell 959, Abell 1351, Abell 1576, Abell 1 722, AbeU 1 995) — gravitational lensing — 
X-rays: galaxies: clusters 
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1. INTRODUCTION 

Rich clusters of galaxies are probably the largest well- 
defined and gravitationally bound objects in the universe, and 
knowledge of their properties is not only interesting in itself, 
but can also set important constrains on models of large-scale 
structure formation. In particular, the cluster mass function has 
caught major interest (e.g., Bahcall & Cen 1993; Gross et al. 
1998). Reliable estimates of cluster masses are also important 
in order to constrain the ratio of baryonic to total mass and to 
determine the density parameter Slo (e.g., White & Frenk 1991; 
Lilje 1992; White, Efstathiou, & Frenk 1993; Carlberg et al. 
1996; Eke, Cole, & Frenk 1996). 

Cluster masses have been measured using different methods, 
having different biases and systematics. The method of clus- 
ter mass measurement with the longest history is the applica- 
tion of the virial theorem to positions and velocities of cluster 
gal axies. The f irst estimates of cluster masses were obtained 



by [Zwicky (1933[ ) and |Smith (1936| ). They found that the total 
mass of galaxies in a cluster only accounted for a small por- 
tion of the total cluster mass. This method works well for low 



redshift clusters, where it can be based on a large number of 
galaxies, but has recently been applied also to galaxies at inter- 
mediate redshifts (e.g., Carlberg et al. 1996) where the number 
of galaxies observed per cluster is usually relatively small. 

Other methods for mass measurement are based on analy- 
sis of gravitational lensing by the cluster and consequent dis- 
tortion of background galaxies and measurements of the X-ray 
temperature of the hot intra-cluster medium. Of these, weak 
gravitational lensing has recently become frequently used for 
estimating the masses of clusters at intermediate and high red- 
shifts (see e.g., recent reviews by Mellier [1999] and Bartel- 
mann & Schneider [2001]). To understand systematic differ- 
ences between the different mass estimators, comparative stud- 
ies are necessary. Dahle et al. (2002; hereafter Paper I) presents 
weak lensing measurements of 39 highly X-ray luminous galax- 
ies at redshift 0.15 < z < 0.35 where the masses and their dis- 
tributions have been determined homogeneously by up-to-date 
methods. In this paper we present velocity dispersion and virial 
mass estimates of five clusters from the sample of Paper I, 
which are compared with the weak lensing results. Our mass 
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determinations ar e based on the method pioneered by L imber 
& Mathews (1960) where the mass estimator is 

My = ^^. (1) 

Here up is the one-dimensional velocity dispersion of the clus- 
ter, while Rh is the projected mean harmonic point-wise sepa- 
ration (projected virial radius) defined by 

N(N- 1) 



Rh = 



(2) 



where N is the number of galaxies in the cluster and R,j is the 
projected distance between any pair of galaxies. In our analysis 
we have used this method with the modifications of C arlberg et 
al. (1996) which are described later in the paper. To be able to 
use equation ([!]) or versions of it for estimating masses, one has 
to both determine the velocity dispersion of a cluster through 
spectroscopy and the angular distribution of cluster galaxies 
from photometric imaging. 

In §2 we present our cluster sample and describe our observ- 
ing and data reduction procedures for the spectroscopic and 
photometric observations of the sample. In §3 we determine 
the velocity dispersions of the clusters, while we in §4 use the 
photometric data to estimate the virial radii of the clusters. The 
virial mass estimates are presented in §5. Finally, in §6, we 
add to our sample a sample of seven clusters for which we have 
obtained weak lensing data and for which there exists veloc- 
ity dispersions in the literature. The velocity dispersions of the 
clusters of this larger sample are compared with our weak lens- 
ing results. 

2. SAMPLE AND OBSERVATIONS 

Our clusters are taken from the weak lensing analysis sam- 
ple of Paper I, and the selection criteria are described in de- 
tail there. That sample consists of 39 highly X-ray luminous 
clusters in the redshift range 0.15 < z < 0.35, taken from the 
samples of Briel & Henry (1993) and Ebeling et al. (1996, 
1998). The six clusters selected for optical spectroscopic analy- 
sis, Abell 914, Abell 959, Abell 1351, Abell 1576, Ab ell 1722, 
and Abell 1995, were all in the Briel & Henry (1993 ) sample. 
Except from Abell 914, they were selected to have X-ray lumi- 
nosities higher than 7 x 10 44 ergs _I in the 0.5-2.5 keV ROSAT 
cluster rest frame band. With the revised redshift presented in 
this paper, also Abell 959 has a somewhat lowe r X-ray luminos- 
ity. The cluster redshifts ( Huchra et al. 1990 ) are in the range 
0.29-0.33, except for Abell 914 which is at z = 0.19. The data 
for Abell 959 were of poor quality because of clouds, and there- 
fore that cluster is omitted from most of our analysis. 

Our spectroscopic observations were made on the three 
nights, 1995 April 30 to M ay 3, with the LDSS-2 Low D isper- 
sion Survey Spectrograph ( Allington-Smith et al. 1994 ) at the 
4.2m William Herschel Telescope (WHT) on La Palma. The 
LDSS-2 was used in a conventional multi-slit mode with slit 
width 1 ."5 and with the medium grism. The slit length was typ- 
ically 15"- 20", which left adequate space for sky line subtrac- 
tion around each galaxy (the half light radii of the galaxies were 
in the range 1"- 1 "5). We used one mask per cluster, with typi- 
cally 30 to 40 slits in the 11' field. The detector was the thinned 
1024 2 pixel Tektronix CCD TEK-1 with an angular scale of 
0"59 pixel -1 giving a dispersion of 5.3Apixel _1 . The FWHM 
seeing was typically 1"1 - l."2 during the observations. The 



spectra normally covered the range 3500 to 8400A, for some 
spectra the range was smaller. 

Target galaxies were chosen based on magnitude and mor- 
phology (color information was not used) from images taken at 
the Nordic Optical Telescope (also at La Palma), coordinates 
were determined from the APM Sky Catalogue 6 , and the slit 
positions of the different masks were then determined in the 
standard way. Each cluster mask was observed with three ex- 
posures, each of 1.8 ks duration, giving a 5.4 ks combined ex- 
posure of each field. Dispersed flat-fields of the twilight sky 
and of a tungsten lamp were taken for each mask. Before each 
cluster observation, an exposure was made of a CuAr wave- 
length calibration lamp with high signal-to-noise lines in the 
range 5400A-7500A. Flux calibration exposures were made of 
a number of spectrophotometric standard stars. 

The data were reduced by the LEXT package (Colless et al. 
1990; Allington-Smith et al. 1994). The data reduction in- 
cluded bias subtraction, correction for distortion effects, flat- 
fielding, combination of exposures, sky subtraction, extraction 
of one-dimensional spectra, and wavelength- and flux calibra- 
tion. The data reduction procedure took into account flexure 
differences between exposures and geometric dis tortions in the 
LDSS-2 optics. The optimal extraction routine of frlorne (1986 ) 
was utilized. The wavelength calibration was based on third- 
order polynomial fits to the CuAr arc lines, and in general pro- 
vided rms residuals less than 0.2A for wavelengths between 
5400A and 8000A, which is the range used to determine the 
redshifts. A typical reduced one-dimensional spectrum of an 
early-type galaxy in one of our clusters is shown in Figure |J 
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FIG. 1 . — Spectrum of an early type galaxy in Abell 1576 with z = 0.302. 

The individual one-dimensional spectra were analyzed with 
the utility crcor which was kindly provided by Dr. W. J. 
Sutherland. Its main estimation method is the cross correla- 



tion technique of Tbnry & Davis (1979), but all spectra are also 
searched for emission lines. In our case, the cross correlation 
was done with a set of 18 template spectra, using galaxy spectra 
from Kennicutt (1992j ) and stellar spectra from Jacoby, Hunter, 
& Christian (1984). The cross correlation routine gives an esti- 
mate of the redshift, the significance of the redshift and its mean 
error, and it also assigns a quality flag ranging from to 4. If 
emission lines are found, a quality flag for the emission line 
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redshift is also assigned based on the number of identified lines 
and their strength. All spectra were also checked manually. The 
cross correlation method was preferred if it gave a quality flag 
equal to or better than the emission line quality flag. In gen- 
eral, redshifts with quality flag 3 or better were found to be 
trustworthy, and were used in the analysis. Spectra with quality 
flag < 2 were discarded, unless they were confirmed by man- 
ual inspection. The data for Abell 959, which were collected 
during partial cloud cover, only yielded five good galaxy red- 
shifts. For the other clusters, the number of good spectra range 
from 13 (Abell 1722) to 24 (Abell 1995). The estimated red- 
shift error was typically 50-100 kms" 1 . One of the objects in 
Abell 1995 was not a galaxy, but a quasar at z = 2.66. Its coor- 
dinates (J2000.0) are a = 14 h 53 m 08.9 s and 6 = 58°03' 1 1 "9. It is 
located about 90" east of the central galaxy of Abell 1995, and a 
simple singular isothermal sphere gravitational lens model with 
a velocity dispersion 1200 kms" 1 (which is consistent with our 
weak gravitational lensing data) implies that the QSO is mag- 
nified by a factor 1.5 by the cluster potential. 

To determine the angular distribution of cluster galaxies, 
photometric imaging of the cluster fields were performed at 
the University of Hawaii 2.24m telescope (UH2.2m) at Mauna 
Kea, Hawaii during three observing runs, 1998 February 19- 
23, 1999, May 14-16, and 2000 March 8-11, using the 8k x 8k 
pixel UH8K mosaic CCD camera in the f/10 Cassegrain focus. 
Further details of these observations and the reduction proce- 
dures of the data are given in Paper I. The UH8K instrument 
has eight 4k x 2k CCDs and has a field of view of 19' by 19', 
enough to fit not only the central parts, but also most of the 
outer regions of the clusters. However, chip # 4 located in the 
NE corner was not used in our analysis, as it had poor cosmetic 
quality. The data were run through a data processing pipeline 
consisting of pre-processing of the original images, transfor- 
mation from chip to celestial coordinates and warping and av- 
eraging to produce the final images. These final UH8K images 
were rebinned into 2k x 2k, giving ."6 pixel" 
the analysis reported in this paper. 
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3. CLUSTER REDSHIFTS AND VELOCITY DISPERSIONS 

In our case, the number of observed galaxies per cluster was 
modest, ranging from 13 to 24 (and for Abell 959 where it was 
not possible to determine a velocity dispersion, five). It is ob- 
vious that estimates of the cluster redshifts and especially ve- 
locity dispersions (and the mass estimator is proportional to the 
square of the velocity dispersion, see eq. ^) must rely on uti- 
lizing the most robust estimators with high resistance to avoid 
problems caused by e.g., deviations from Gaussian velocity dis- 
tribution and contamination from foreground and background 
galaxies. In our analysis, we used the optimal estimators for 
contaminated galaxy samples of this size (Beers, Flynn, & Geb- 
hardt 1990), the biweight estimators of location and scale, as 
estimators for the cluster redshifts and velocity dispersions, 
respectively. Confidence intervals were calculat ed using the 
bias corrected and accelerated bootstrap method (Efron 1987) 



with 10,000 bootstrap replications. These robust estimators and 
methods for determining confidence intervals have been well 
tested on redshift samples like ours (except from Abell 959, see 
e.g., Beers et al. 1990; Girardi et al. 1993; Borgani et al. 1999). 
For the biweight estimator of location, we chose a value for 
the tuning constant (c = 6.0) that gave the estimator good ro- 
bustness and high efficiency. For example, in the case of a sam- 
ple approximating a Gaussian distribution, the rejection point is 



cMAD = 4 standard deviations from the median, where MAD 
(= 0.6745 in the Gaussian case) is the median absolute devia- 
tion from the sample median (for details, see e.g., Mosteller & 
Tukey 1977). For the biweight estimator of scale, we chose a 
tuning constant c = 9.0 such that in the Gaussian case, galaxies 
more that 6 standard deviations away from the median would 
be rejected. Our statistical analysis was performed using the 
ROSTAT program kindly provided by Dr. T. C. Beers. A single 
iteration was performed for the biweight measurements. 

Our redshift and velocity dispersion estimates are presented 
in Table [j]. The analysis uses only galaxies with well measured 
redshifts which have not been rejected as foreground or back- 
ground objects. The second column in Table |I| N, shows the 
number of galaxy redshifts used in each cluster. The third and 
fourth columns show the derived cluster redshift and the derived 
one-dimensional velocity dispersion together with 68% confi- 
dence limits. Given the modest number of spectroscopically 
measured galaxies in each cluster, we performed an additional 
estimate of the errors associated with small-number statistics to 
check if they were compatible with the confidence intervals in 
Table |l| This was done by repeatedly drawing random samples 
of galaxies (with N =10 and N = 20) from a much larger redshift 
data set (N ~ 200; see Yee et al. 1996) for Abell 2390 (z = 0.23), 
which has a velocity dispersion of <jp = 1093 ± 61 kms" 1 (Carl- 
berg et al. 1996). The variance in the values for the redshift 
and velocity dispersion calculated from 100 such subsamples 
is similar to, or slightly (< 30%) larger than, the confidence 
intervals in Table |l[ 
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FIG. 2. — Color-magnitude diagram for Abell 1351. The positions of the 
galaxies with measured redshifts and colors are plotted with open circles. With 
a few exceptions, the selected galaxies are all within the main cluster sequence. 

For most of our six clusters, the determined redshifts are in 
relati vely good agreement with previous results (H uchra et al. 
1990), but ours are based on a much larger number of galaxy 
redshifts. However, for Abell 959 where we only base our 
cluster redshift on four g alaxies, our redshift is significantly 
lower than the redshift of |Huchra et al. (199"0 ). Three of the 
four galaxies that we based our redshift estimate on have V -I 
colors putting them in the middle of the red sequence of the 
col or-magnitude diagr am of the field. Both the galaxies used 
by |Huchra et al. (1990 ) were substantially bluer than the cluster 
red sequence. We therefore believe that our cluster redshift for 
Abell 959 is correct. The 68% confidence limits for the redshift 
of this cluster could not be reliably estimated using the boot- 
strap method we used for the other clusters, and the confidence 
limits shown in Table [j] are instead based on the variance in the 
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Table 1 

Cluster redshifts and velocity dispersions 





N 


z 


up (kms L ) 


z (H90) a 


N (H90) b 


Abell914 


11 


1934+0-0017 


1140^ 


0.1941 


3 


Abell 959 


4 C 


/_o. 003 4 




0.3530 


2 


Abell 1351 


17 


397Q+0.0014 


1680+g , 


0.3220 


3 


Abell 1576 


14 


2986+ 00011 


1040!j™ 


0.3020 


1 


Abell 1722 


9 


^964+0 0013 

u.jzof_ 0()017 


970ig° 


0.3275 


1 


Abell 1995 


20 


3207+ 00010 


113<£i?° 


0.3180 


2 



"Literature redshift ( Huchra et al. 199C ). 



Huchra et al. (1990). 



'Number of individual galaxy redshifts measured by I 
' With only four cluster galaxies with measured redshifts in the Abell 959 field, a velocity dispersion could not be estimated. The confidence interval for the cluster 
redshift was estimated from the variance in the results obtained by randomly drawing redshift values from a much larger data set for a similar cluster (see the text for 

details). 



results obtained by drawing 4 galaxies at random from the Yee 
et al. (1996) data set for Abell 2390. 

Struble & Rood (1999) list a redshift of z = 0.279 for 
Abell 1576, significantly lower than the value we derive from 
our data and the value measured by Huchra et al. (1990 ). The 
lower redshift value (which is also the one given by NED) orig- 
inates from a study by Leir & van den Bergh (1977) and is an 
estimate based on a combination of cluster galaxy photometry 
and cluster richness. We note that Abell 1351, where the mea- 
surement is based on 17 redshifts, has an abnormally high ve- 
locity dispersion. This is partially confirmed by the weak lens- 
ing analysis (see §6). Patel et al. (2000) have recently published 
redshift and velocity dispersion measurements of Abell 1995, 
based on the redshifts of 15 galaxies, six of which are com- 
mon with our data set. Their result of z = 0.322 ± 0.001 and 
(Tp = 1282^}20kms _I agrees quite well with our result. 

4. ANGULAR DISTRIBUTION OF CLUSTER GALAXIES 

As described in §1, the virial mass estimator (eq. [I]) tradi- 
tionally uses the projected mean harmonic point-wise separa- 
tion defined in equation (Q). However, this estimator is sensi- 
tive to close pairs and is somewhat noisy (see e.g., Carlberg et 
al. 1996). It will also systematically underestimate the radius 
for a rectangular aperture wh en the entire cluster is not covered 
(Bahcall & Tremaine 1996). Carlberg et al. (1996 ) proposed to 
instead use a related radius, the ring-wise projected harmonic 
mean radius, defined by 

N(N— 1) 



Rh 



J_ r 2 * de 

<J 27r J0 y/Rj+Ri+UtiRjCose 



T — f 
N(N- 1) 



(3) 



£,-< 



-,K(k u ) 



i<j ir{Rj+Rj) 

where Rj and Rj are the distances of galaxies i and j from 
the cluster center, kfj = 4R i R j /(R i +Rj) 2 and K(k) is the com- 
plete elliptic integral of the first kind in Legendre's notation. 
This radius treats one of the particles in the pairwise poten- 
tial |R, -R, | _1 as having its mass distributed like a ring around 
some given cluster center and is less noisy and handles non- 
circular apertures better than the standard Rh- It will, however, 
overestimate the true projected virial radius if the cluster is sig- 
nificantly flattened or extensive sub-clustering is present. In our 
analysis we therefore use both methods to check the influence 



of the choice of radius in the virial mass estimator on the un- 
certainty of the resulting virial mass. 

The radii R/, and Rh should in principle be computed for all 
galaxies residing in the virialized regions of each cluster. In 
our case, we have only measured redshifts of a few galaxies 
in each field, and redshifts can therefore not be used to select 
likely cluster galaxies over which these sums are to be taken. 
Instead we use the UH8K photometry described in §2 to se- 
lect the cluster galaxies. Our procedures for object detection 
and classification follow Kaiser, Squires, & Broadhurst (1995), 
and are described in detail in Paper I. Briefly, we first ran a 
peak finding algorithm with a 4er detection limit. Then we 
performed aperture photometry on the detected objects and re- 
moved stellar objects brighter than / = 21'" by applying cuts in 
the size-magnitude diagram. After having removed stars and 
very faint objects (the completeness limit is / ~ 25.0'"), the V- 
and /-catalogs were combined. In the final catalog we kept ob- 
jects that were detected in both V and / with positions that dif- 
fered by less than l."2 (2 pixels). The peak finding algorithm 
of Kaiser et al. (1995) is optimised for detecting small, faint 
galaxies, and may occasionally have problems resolving closely 
grouped galaxies in cluster cores. To test whether our choice of 
galaxy detection algorithm had a significant influence on our 
results, we also calculated Rh values from galaxy catalogs gen- 
erated using the SExtractor object finding algorithm (Bertin & 
Arnouts 1996). The resulting variations in the Rh values were 
small and within the quoted uncertainty intervals. In addition, 
the performance of the Kaiser et al. (1995) detection algorithm 
was tested by manually examining images of cluster fields with 
the catalog object positions superposed as small circles. The 
positions indicated by the Kaiser et al. algorithm were indeed 
almost always centered on the brightest galaxy within the circle 
(for an example, see Fig. |]), and we could not manually find 
obvious cluster members not found by the algorithm. 

In a color-magnitude diagram (see Fig. | for Abell 1351), the 
cluster galaxy red sequence consisting of early type galaxies is 
clearly visible (see also e.g., Olsen [2000]). We exploited this 
fact in our procedures for discriminating cluster members from 
galaxies possibly belonging to background or foreground struc- 
tures. As described in detail in Paper I, we generated catalogs 
of "cluster" galaxies by making appropriate cuts in each clus- 
ter color-magnitude diagram to include early type galaxies at 
the cluster redshift. We also tried to use galaxy samples that 



Cluster virial masses 



5 




FIG. 3. — Image of Abell 1351 overlaid positions of cluster galaxies identified for the color-magnitude diagram (circles). In this figure, north is up, east is left and 
the scale on the axes is in pixels of size 0"6. 
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included all galaxies in the fields down to a certain magnitude 
limit without any color selection. This did not produce sig- 
nificant changes in the results, except for Abell 1351, which 
displays clear evidence of a foreground structure (the color- 
magnitude diagram for this cluster clearly shows an additional 
red sequence at V— I ~ 1.2). To assess the uncertainty in the 
analysis caused by some arbitrariness in the exact positions of 
the color cuts, we varied the limits of the cuts to make an en- 
semble of galaxy catalogs for each cluster on which the further 
analysis was performed in parallel. To determine Rh, the po- 
sition of the cluster center is needed. This was found in two 
steps, first by applying a peak finding algorithm on a smoothed 
map of the distribution of galaxies defined to be cluster mem- 
bers by the color criteria described above. As described later 
in this chapter, we chose three different sets of galaxies within 
three different radii from this center to perform our analysis. 
For each of these sets we recomputed the center which we used 
for determining R/,. The center was set to be the biweight esti- 
mate (Beers et al. 1990) of the mean of the x and y coordinates 
of all galaxies used in the estimate. We tested for the effect of 
varying the center position, and found that varying this by less 
than 50 pixels did not significantly affect the results. 

The radius entering the virial mass estimator (eq. [[]), and 
hence the mass estimate itself, can be biased if we include 
galaxies outside the virialized part of the cluster, or vice versa, 
if we exclude galaxies outside some radius within the virialized 
part of the cluster. The projected density distribution of galax- 
ies outside the cluster core is well known to fall approximately 
oc R' 1 both in the virialized region and in the surrounding in- 
falling region (reflected in the cluster-galaxy cross-correlation 
function [e.g., Lilje & Efstathiou 1988]). Thus we would not 
expect to see a visible border between the virialized region and 
the infalling region. This is confirmed by our data. However, it 
is normally assumed (see e.g., Carlberg et al. 1996) that the viri- 
alized region encompasses the region around the cluster center 
where the average density is about 200 times the critical density 
of the universe. This is based on the fact that in the Einstein de 
Sitter universe, a purely spherical top-hat perturbation will viri- 
alize at a density 178 times the background density at the time 
of virialization. In other universe models, the value is some- 
what different, and in open models all clusters virialize at an 
epoch much earlier than the present epoch. However, the fig- 
ure 200 is sufficient as an approximate estimate (see e.g., Lilje 
1992). For a typical rich cluster this is the region within a ra- 
dius quite close to the Abell radius (1 .5/r'Mpc). Our 19' by 19' 
field is about two Abell radii across at the typical cluster red- 
shifts, so our field is approximately of the right size to minimize 
this possible bias. In our analysis, we compute Rh and R/, both 
using all "cluster" galaxies in the field, and for only the galaxies 
within radii of 6 respectively 8 arcminutes from the cluster cen- 
ter. The scatter between these three different estimates of the 
radii estimated from the same cluster was normally of the order 
of 10%. There was a systematic effect, in that the estimated 
radii to some extent scaled with the size of the cutoff. We use 
the differences in these results to determine the uncertainty in 
our final virial mass determinations. As seen in Figure ^ where 
the positions of all "cluster galaxies" in a typical catalog for 
Abell 1351 are shown as circles, most of the galaxies defined 
as "cluster" galaxies in each field do indeed reside within the 
inner of these radii. Before doing our analysis we also had to 
mask certain areas in the fields and reject galaxies in those ar- 
eas from our analysis. Those areas included the UH8K CCD 



# 4 and the areas surrounding bright stars. For Abell 914, the 
cluster was not positioned in the center of the field but in the 
NW quadrant, to be able to simultaneously image another clus- 
ter in the SE quadrant. For Abell 914 we therefore only used 
galaxies in the NW 15' x 15' part of the mosaic in our analysis. 

For each cluster we constructed three ensembles of "cluster" 
galaxy catalogs, using the whole image, the central 8' field and 
the central 6' field. Each ensemble consisted of about 5 catalogs 
where the color and magnitude selection criteria for "cluster" 
galaxies had been varied within reasonable bounds. For each 
catalog in all these ensembles we computed Rh (eq. ||) and Rj, 
(eq. |J). For all the ensembles, Rh differed from Rh by less than 
5%, and we chose to only use Rh in our further analysis. We 
express R/, in length units assuming an Einstein de Sitter cos- 
mology (Slo = 1, A = 0). If we instead had assumed a cosmology 
with f^o = 0.3 and A = 0.7, the radii would have become 13% 
larger at z = 0.3 and 9% larger at z = 0.2. Compared to our other 
uncertainties, this error is relatively small. Finally, we express 
the radius used in our further virial mass analysis as the depro- 
jected virial radius, 



7T 

rv = ^R h . 



(4) 



Dr. T. C. Beers' ROSTAT program (described in §3) was then 
applied to the combined ensemble of estimates of ry for each 
cluster to find the best estimate of its ry (by the biweight esti- 
mator of location), and the 68% confidence limits on this esti- 
mate (by the bias corrected and accelerated bootstrap method). 
With the systematic effect mentioned above, i.e., a correlation 
between the estimated R/, and the radius within which galaxies 
have been picked for analysis, it is clear that this method has its 
deficiencies. It does not necessarily give an unbiased estimate 
of the expectancy value of ry, and the estimated confidence lim- 
its cannot be treated as formal 68% confidence limits. However, 
with the rather large uncertainties in these estimates, the quoted 
confidence limits give a fair estimate of the uncertainty in the 
results, which are of order 10%. Our results are presented in the 
second column of Table [l| Since the small number of redshifts 
for Abell 959 prevented us from finding its velocity dispersion, 
we also omitted it from the galaxy distribution analysis. 

Table 2 

Cluster Virial Radii and Masses 





ty (h 'Mpc) 


M v (lO'^/r'Mo) 


Abell 914 


1 23 +,UJ 


1.1 ±0.6 


Abell 1351 
Abell 1576 
Abell 1722 
Abell 1995 


L6 8 |§ 
1 52+°- 10 

1.80 „ , g 


3.3 ± 1.4 

4.0 ±1.4 
1.2±0.6 

1.1 ±0.4 



5. CLUSTER VIRIAL MASSES 

The virial masses of the clusters, given by equation (jl|), are 
shown in the third column of Table 0. The given statistical er- 
rors are based on half the length of the 68% confidence intervals 
of the cluster velocity dispersions and virial radii, and standard 
propagation of errors. 

The virial mass given by equation ([j]) and Table || is the mass 
of the virialized region of the cluster, supposed to be within the 
virial radius ry. Since the density profile in the outer parts of 
the virialized region, and at least the inner part of the infalling 
region, in a cluster is close to being oc r~ a where 2 < a < 3 
(e.g., Adami et al. 1998; Carlberg et al. 1996; Lilje & Efstathiou 
1988), this mass can be scaled to other radii assuming that the 
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mass inside a given radius grows approximately proportional to 
the radius. 

As pointed out in §4, there is no clear feature seen in clusters 
at the transition from the virialized region to the surrounding in- 
falling region. Also, the idea of a virialized region bounded in 
space by a spherical border, surrounded by the non-virialized 
infalling region, is clearly an oversimplification. There might 
well be recently infallen galaxies or subclusters within the virial 
radius which have not been virialized, and these will bias the 
mass estimate. Also, we know that real clusters may be far from 
spherical. Even the quite large statistical errors given in Table ||, 
should therefore be taken with a pinch of salt, because of the un- 
certainty in what the derived virialized mass really measures. 

Even for our idealized model, there are also several other 
possible sources for systematic errors. If our measurement only 
uses galaxies inside a bounding radius much smaller than the 
virial radius, our mass estimate will be an overestimate of the 
mass within r/, be cause of the neglect of the surface term in 
the virial theorem (Carlberg et al. 1996). This effect is at most 
50%. Our virial radius estimates of §4 were derived assuming 
an Einstein-de Sitter cosmology. As pointed out in that chapter, 
this overestimates the virial radius, and hence the virial mass, 
by about 10% if the real universe is vacuum energy dominated 
with n ~ 0.3. 

6. COMPARISON WITH WEAK LENSING MEASUREMENTS 

In Paper I we derived weak lensing mass estimates for a 
sample of 39 galaxy clusters including the clusters studied in 
this paper. Seven additional clusters from the Paper I sam- 
ple have spectroscopically determined velocity dispersion mea- 
surements available in the literature. Here, we compare these 
measurements to the velocity dispersion estimates ctwl derived 
in Paper I by fitting a singular isothermal sphere (SIS) profile to 
the observed tangential weak lensing distortion as a function of 
radius. We note, however, that three of the clusters (Abell 959, 
Abell 1351 and Abell 1576) were poorly fit by a SIS profile. 

In Table | we summarize the available data for the clusters 
discussed in this paper, and the available data for the clusters 
with published velocity dispersion measurements for which we 
have weak lensing data. 

The quantity odm (= &wl listed in Paper I) in the fourth col- 
umn of Table is derived from the weak lensing results of Paper 
I. A subsample of 1 1 clusters have spectroscopically measured 
velocity dispersions ap, which are plotted against the dark mat- 
ter velocity dispersions ctdm in Figure [5J. 

There is generally a reasonably good agreement between ap 
and itdm, with a tendency for the former to be smaller than 
the latter by a barely significant amount, the median value of 
cf/ctdm being 0.94. We also find (a P /a DM ) = 1.02 ±0.13 
when using all twelve clusters and (o> / odm) = 0.97 ± 0.08 af- 
ter excluding two outliers mentioned below. If we only con- 
sider the clusters for which we have measured <jp in this paper 
(but excluding Abell 914, for reasons given below), we find 
(o>/Vdm) = 0.98 ±0.12. 

One cluster (Abell 914, at z = 0.194) show a<T P a factor — 2 
larger than (Tdm, with a very strong (> 2a) discrepancy between 
the two velocity dispersion values. However, there is a second 
cluster (Abell 922, at z = 0.189) at a projected distance of only 
11 arcminutes on the sky. Since the weak lensing measure- 
ments indicate that Abell 922 is a more massive cluster than 
Abell 914 (its estimated odm is SlO^J^kms" 1 ), the observed 
redshifts, which were measured in a 10 x 10' field centered on 
Abell 914, may have been significantly contaminated by galax- 



ies along the line of sight belonging to the other cluster. 

A second cluster (Abell 697, at z = 0.282) shows an al- 
most equally strong discrepancy in the opposite direction, with 
f£)M ~ 2 x a p. We note however, that the spectroscopically 
measured galaxy velocity dispersion is based on a rather small 
number of cluster galaxies (N = 9; see Metzger & Ma 2000). 
o 
o 
o 




500 



1000 



1500 



2000 



a nM [km s ] 



FIG. 4. — The spectroscopically measured velocity dispersion ap vs. the 
dark matter velocity dispersion <tdm- The error bars shown are at Icr. The 
dotted line indicates equality between the values. Filled circles represent the 
clusters for which we present ap measurements in this paper. The other points 
have ap values drawn from the literature. 



7. DISCUSSION 

Our in general good agreement between ap determined from 
dynamical estimates and ap,M determined from gravitational 
lensing estimates contrast other results. 

Smail et al. (1997) presented a comparison of weak lensing 
and spectroscopic mass measurements for a sample of 9 clus- 
ters at 0.17 < z < 0.55 which were imaged with the WFPC2 
instrument aboard the Hubble Space Telescope (HST). They 
compared the average tangential shear observed within an an- 
nulus around the optical cluster center (60 < r < 200/r'kpc) 
to the expected shear from SIS-type mass profile with the ob- 
served galaxy velocity dispersion. The shear signal was found 
to be significantly lower than predicted from the spectroscopic 
measurements for all the clusters, and Smail et al. (1997) sug- 
gested that this was caused by systematic overestimates (by typ- 
ically ~ 50%) of the cluster velocity dispersion values given in 
their paper. As noted in §| we see no such discrepancy in our 
data. Hoekstra et al. (2001) argue that this discrepancy is in- 
stead a result of the small field size of the single WFPC2 point- 
ings used by Smail et al. (1997) to measure the gravitational 
shear. If there is significant substructure in the central regions 
of the clusters, or if the density profiles of the clusters are sig- 
nificantly shallower than a SIS-type model at small radii, the 
shear measurements are likely to underestimate the true cluster 
mass. 

Brainerd et al. (1999) have used ray-tracing through high- 
resolution TY-body simulations of massive (My ~ lO^/r'M©) 
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Table 3 
Cluster data. 



Designation 


z 


Op 


Ker. 




r 

Lx 

^W. 1 — JvC V J 






(Vm <T ls l 


Up 




10 44 prac -1 


AK^II 1 1 ^ 
/\Deii 1 1 j 


1 Q7 


1U/4_ 12] 




iiju_ 270 


1 /I CQa 


/\Dcll jZU 




700 It /O 




1 rwn 4- 1 nn 
llOU it 1UU 


l^f . JZ 


Audi 003 


U. to/ 


091+233 




1010_ I70 




A Kail 
ADell 0y 1 




Q/l 1 _i_ 90/; 




173(J -200 


in icb 


Ahfll Q14 




1 1 4.1 + 153 






s 00 b 

J -v/W 


Abell 959 c 


0.286 






990!|°° 


14.3 b 


Abell 963 


0.206 


1350!?°° 


(4) 


1070!{| 


8.72 b 


Abell 1351 


0.328 


1680 + ^o 




1410^0 


8.31 b 


Abell 1576 


0.299 


1041!!| 




1060 ±90 


11.52 b 


Abell 1722 


0.326 


966l?g 




1160 ±140 


9.70 b 


Abell 1995 d 


0.321 


1 1 ?6 +I51 




1240 ±80 


13.42 b 


Abell 2104 


0.153 


1200 ±200 


(5) 


1390±180 


7.89 e 


Zwicky 7160 


0.258 


1133 ±140 


(2) 


1230±130 


13.73 b 



Note. — The values are given for an Einstein-de Sitter Universe. See the text for a brief discussion of possible additional systematic uncertainties that are not 
included in the error intervals given here. 

References. — (1) Girardi & Mezzetti 2001, (2) Carlberg et al. 1996, (3) Metzger & Ma 2000, (4) Lavery & Henry 1998, (5) Liang et al. 2000 

"X-ray luminosity measured by Ebeling et al. (1998). 
*X-ray luminosity measured by Bohringer et al. (2000). 
' The values have been corrected for the revised redshift for Abell 959 presented in this paper. 
rf Patel et al. (2000) measure a P = 1282+[g for Abell 1995. 
c X-ray luminosity measured by Ebeling et al. (1996). 



clusters to show that the average tangential shear, when mea- 
sured within an annulus with an outer radius Rm/kx significantly 
smaller than the virial radius (and assuming a SIS-type mass 
profile), will tend to underestimate the true cluster mass. For 
Rmax — Rv, they find that the SIS-model fit yields quite accu- 
rate mass values (typically 5% - 10% underestimates; see their 
Figure 9) even if the true mass density profile of the clusters is 
more accurately represented by the model proposed by Navarro, 
Frenk, & White (1995). We note that for most of the clusters 
listed in Table ||, we have estimated a DM based on tangential 
shear measurements out to radii ~ Rv (for most of our clusters, 
we have used Rmax = 550" = 1 .5/T 1 Mpc at z = 0.3; see Paper I 
for details), so we expect our values to be essentially unbiased. 

Girardi & Mezzetti (2001) have reanalyzed the available ra- 
dial velocity measurements in the clusters studied by Smail et 
al. (1997) and derive velocity dispersion values that are on av- 
erage about 20% lower than those given by Smail et al. (1997). 
Thus, it appears that the discrepancy between spectroscopic 
and lensing mass estimates seen by Smail et al. (1997) is most 
probably caused by a simultaneous overestimation of dynam- 
ical masses and underestimation of lensing masses, of almost 
equal magnitude. Girardi & Mezzetti (2001) also compare their 
virial cluster mass estimates with weak lensing mass estimates 
drawn from the literature for 18 clusters (for 7 of these the lens- 
ing mass was measured from single HST WFPC2 pointings), 
and find a median value M spec /Mi em = 1.30 and a 90% confi- 
dence interval 0.63 < M spec /Mi em < 2.13. This is consistent 
with our results in §|| within the errors, but the Mi ens values 
they use are probably subject to a larger bias caused by the (on 



average) significantly smaller fields used to measured the shear. 

In their study, Brainerd et al. (1999) did not consider the pos- 
sible effect of structures outside their clusters on the lensing 
mass estimates. Galaxy clusters are believed to be connected 
by large filaments of dark matter, diffuse gas, and galaxies, and 
mass within such filaments will contribute significantly to the 
lensing mass when seen in projection along the line of sight to 
a cluster, even if it is physically separated from the cluster by 
several Mpc (Metzler, White, & Loken 2001). Still, even in 
this scenario, odm may not necessarily be significantly larger 
than up, since the galaxy velocity dispersion would also be bi- 
ased upwards by galaxies outside the virial radius falling into 
the main cluster along the filaments. 

The modest fraction (~ 30%) of massive clusters that have 
significant substructure in their dark matter distribution provide 
evidence that the majority of massive clusters at z = 0.2-0.3 
are close to dynamical equilibrium (Dahle et al. 2002a). The 
good agreement that we find between virial masses and lensing 
masses for most of our clusters is consistent with this picture, 
although significant amounts of substructure could probably be 
present in the clusters without changing the ap values beyond 
the error bounds of our measurements. To reliably assess the 
dynamical state of the clusters from spectroscopic data alone, 
it would probably be necessary to observe galaxy samples an 
order of magnitude larger than those we present in this paper. 
In addition, X-ray observations, particularly using Chandra and 
XMM, would be an efficient aid to measure substructure and 
determine the dynamical state of the clusters. We plan to inves- 
tigate the relation between weak lensing and X-ray properties 
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for a larger sample of clusters in a future paper (H. Dahle et al. 
2002b, in preparation). 
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